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31 P and 75 As NMR evidence for a residual density of states at zero energy 
in superconducting BaFe 2 (As .67Po.33)2 
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31 P and 75 As NMR measurements were performed in superconducting BaFe2(Aso.67Po.33)2 with 
T c = 30 K. The nuclear-spin-lattice relaxation rate Tj and the Knight shift in the normal state 
indicate the development of antiferromagnetic fluctuations, and in the superconducting (SC) 
state decreases without a coherence peak just below T c , as observed in (Bai_ a; K a; )Fe2As2. In contrast 
to other iron arsenide superconductors, the oc T behavior is observed below 4K, indicating the 
presence of a residual density of states at zero energy. Our results suggest that strikingly different 
SC gaps appear in BaFe2(Asi_a;Pa:)2 despite a comparable T c value, an analogous phase diagram, 
and similar Fermi surfaces to (Bai„ a; K :E )Fe2As2. 
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The symmetry of the superconducting (SC) order pa- 
rameter is one of the most important issues in iron 
pnictide superconductors^ NMR studies on i?FeAsO 
("1111", R is a rare-earth element) 2 ^— and ^4Fe2As2 su- 
perconductors ("122", A is an alkaline-earth element 
)11±11 have consistently found the absence of a coherence 
peak, but T-depcndences of the spin-lattice relaxation 
rate Tf 1 below T c remain controversial. Although these 
various behaviors of below T c may reflect that Tj -1 
in iron pnictidcs is sensitive to doping levels and/or sam- 
ple quality, these NMR results 2 - - — show no evidence of a 
residual density of states (DOS) at zero energy in the SC 
stated This is in sharp contrast to the case of high-T c 
cupratcs, where dilute impurities and crystalline defects 
easily induce a residual DOS at the Fermi level. 

Recently, the penetration depth study on the iron 
phosphide LaFePO (T c ~ 6 K) suggested nodal super- 
conductivity^ which is inconsistent with the prevailing 
theory of fully-gapped s± states*^— This different SC 
state is shown to be explained by changes in pnictogen 
height acting as a switch between high-T c nodeless and 
low-T c nodal pairing^ In this regard, the isovalent P- 
substitutcd BaFc2(Asi_ 2; P a ;)2 where P substitution is ex- 
pected to reduce the pnictogen height is interesting in 
order to investigate the nature of the SC gap function. 
Isovalent P substitution for As can tune the ground state 
of BaFe2As2 from antiferromagnetism to superconduc- 
tivity with a comparable maximum T c value of 30 K, 
similar Fermi surfaces, and a similar phase diagram to 
(Bai_ K K ;c )Fe2As2 (see Fig. 1) 22*21 It is believed that 
P substitution in BaFe2(Asi_ a; P a; )2 does not introduce 
charge carriers, unlike the hole-doped (Bai_ :c K a ;)Fe2As2 
and electron-doped Ba(Fei_ s Co x )2As2. Here, we present 
NMR evidence for a residual DOS at zero energy in 
BaFe2(Aso.67Po.33)2, which is fundamentally different 
from other iron arsenide superconductors. Our findings 
suggest that different SC pairing states are realized in 
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FIG. 1: (Color online) ac susceptibility of BaFe2(Aso.67Po. 33)2 
obtained with an NMR coil at ~ 17 MHz, indicating 
T C (H = 0) = 30 K. Inset: schematic phase diagram 
of BaFe2(Asi_ a: P :E )2. To (Tafm) represents the structural 
(AFM) ordering temperature. The arrow denotes the P sub- 
stitution x of the present sample with maximum T c . 



iron pnictides despite comparable high-T c values, analo- 
gous phase diagrams, and similar Fermi surfaces. 

Polycrystallinc samples of BaFe2(Asn.67Po.33)2, con- 
sisting of single crystals with average dimensions of 100 
x 100 x 50 fun , were prepared from mixtures of FcAs, 
Fe, P, and Ba^ The P substitution x was determined 
by an energy dispersive x-ray. ac-susceptibility measure- 
ments with an NMR coil show that T C (H = 0) is 30.0 K 
with a narrow transition width, indicating an excellent 
sample quality (see Fig. 1). 

The inset in Fig. 2 shows 31 P NMR spectra obtained 
by sweeping frequency at a fixed magnetic field (4.12 T). 
The Knight shift K was determined with respect to a 
resonance line in H3PO4 (K = in Fig. 2). K is T- 
independent (K = 0.131 ± 0.001%) in the normal state 
and decreases in the SC state. From the magnetic field 
dependence of K, we find that this decrease is dominated 
by the SC diamagnctic shift in the mixed state. The 



2 



T (4.12 T) 



<^CO 



\ A 





• 




A 




• 


• 


A 


• 






A 




6.00 T 


• 


4.12 T 




2.10 T 




71.0 71.2 71.4 
f(MHz) 



T(K) 

FIG. 2: (Color online) 31 P Knight shift K for different mag- 
netic fields. Inset shows 31 P NMR spectra in 4.12 T. The 
spectra have been offset vertically. 
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FIG. 3: (Color online) 31 P (right axis) and 75 As (left axis) 
(TiT)" 1 for BaFe 2 (Aso. 6 7Po.33)2, along with 75 As (TiT) -1 in 
optimally doped iron arsenide superconductors.— 1 ^ 1 10 ' 12 i 14 i 23 

large SC diamagnetic shift prevents us from determining 
variations in the spin susceptibility in the SC state. The 
31 P NMR spectra show no evidence for antiferromagnetic 
(AFM) order at any temperature.— 

The 31 P ( 75 As) nuclear spin-lattice relaxation rate 
T^ 1 was determined by fitting the time dependence of 
the spin-echo intensity to a theoretical I = 1/2 (3/2) 
curve after a saturation pulse. Good fits are obtained 
throughout the temperature range 100 mK - 270 K (see 
Fig. @|, indicating a homogeneous electronic state in 
BaFe2(Aso.67Po.33)2- In contrast to the T-indcpcndcnt 
Knight shift, (TiT) -1 for 31 P and 75 As increases signifi- 
cantly on cooling down to T c (see Fig. [3]). Since (TiT) -1 
is related to the sum of q-averaged dynamical spin sus- 
ceptibility x(q) a t low energy and the Knight shift is pro- 
portional to x(q = 0), this (TiT) -1 result clearly shows 
the development of low-lying spin fluctuations at finite 
wave- vectors Q away from zero. 

Next, we compare spin dynamics in 
BaFe 2 (As .67Po.33)2 with those in (Bao. 6 Ko.4)Fe 2 As2, 



Ba(Fco.9 2 Coo.o8)2As 2 , LaFeAs(O .89F .ii), and 
LaFeAsOo.6 superconductors, which are all near 
optimum doping. As in BaFe 2 (Aso.67Po.33)2, the devel- 
opment of low-energy AFM fluctuations is inferred from 
an increase in (TiT) -1 in (Bai_xKa;)Fe2As2 (Refs. 10- 
13)^1 and Ba(Fei_ x Co x ) 2 As 2; ±^ suggesting the im- 
portance of AFM fluctuations for superconductivity (see 
Fig. 3). In contrast, low-lying AFM fluctuations are sup- 
pressed above T c in optimally doped LaFcAs(Oo.89Fo.n) 
(Refs. 2 and 3)2^ and LaFeAsO .6r^ implying different 
nature of low- lying AFM fluctuations in 122 and 1111 
superconductors. Interestingly, the T-dependences of 
(TiT) -1 differ even among the Bal22 superconductors; 
(TiT) -1 in (Bao.6Ko.4)Fc 2 As 2 exhibits the most signif- 
icant T-dcpendence, there is only a slight enhancement 
below about 50 K in Ba(Feo.9 2 Coo.os)2As 2 , and an 
intermediate T-dependence of (TiT) -1 is observed in 
BaFe 2 (Aso.67Po.33)2- These variations indicate that 
low-lying AFM fluctuations are sensitive to composition, 
and this is attributable to different nesting conditions in 
their Fermi surfacesi 22 ! 25 Although these materials have 
similar Fermi surfaces with three hole sheets and two 
electron sheets, band calculations suggest that P and 
K substitutions would modify their Fermi surfaces in 
different way a 22 ' 25 ; P substitution makes the hole sheets 
having stronger c-axis warping, and thus nesting is 
only possible in restricted regions along the c direction, 
while K-doping makes the hole sheets larger and the 
electron sheets smaller, resulting in only one electron- 
and hole-sheet among five contributing to the nesting. 

Figure 3 highlights systematic variations in (TiT) -1 
from hole- to electron-doped arsenides. The large values 
of (TiT)- 1 in the hole-doped (Ba . 6 K .4)Fe 2 As 2 at high 
temperatures most likely originate from the enhanced 
DOS at the Fermi level inferred from band calculations, 25 
In contrast, in the electron-doped Ba(Fco.9 2 Coo.os)2As 2 , 
LaFeAs(0 .89Fo.n), and LaFeAsO .6- (TiT) -1 is smaller 
and decreases on cooling at high temperatures. This is 
attributable to a hole Fermi surface with a high DOS 
just below the Fermi levels Thus, the variations in 
T-dependences of (TiT) -1 suggest that spin dynamics 
probed by NMR arc strongly related to the nature of the 
Fermi surfaces. 

In the SC state of BaFe 2 (Aso.67Po.33)2, we observed a 
marked Tj -1 behavior distinctly different from other iron 
arsenides. Figure |3] (c) shows that 31 P T^ 1 decreases 
sharply just below T c and there is no evidence for a co- 
herence peak. We find Tf 1 oc T below ~ 4 K down to 
100 mK, showing the existence of a residual DOS at zero 
energy. Below about 10 K (« 0.3T C ), T x -1 depends on 
the magnetic field [sec Fig. 4 (c)]. By extrapolating to 
zero field, we still find T]T = const behavior below ~ 4 
K, excluding the possibility that the TiT = constant be- 
havior is attributable to applied fields. The existence of a 
residual DOS in H — > inferred from NMR is consistent 
with thermal conductivity results on BaFc 2 (Aso.67Po.33)2 
that reveal a finite thermal conductivity divided by tem- 
perature kq/T at T — > under zero fields Our findings 
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FIG. 4: (Color online) (a) and (b) show 31 P nuclear-spin re- 
laxation curves for 175 K and 1.5 K in 4.12 T, exhibiting clear 
single exponential behavior, (c) 31 P T^ 1 under 2.1, 4.12, and 
6.0 T. Inset in panel (c) shows 31 P Tf 1 vs magnetic field at 
2 K; a linear fit is used to extrapolate T^ 1 for H — > T. 
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FIG. 5: (Color online) Normalized 31 P T^ 1 for 
BaFe 2 (As .67Po.33)2 and (Ba .6K .4)Fe 2 As2^ The solid 
curve shows a calculation for a line-node gap function 
A(4>) = A sin (2^) with 2A = 6k B T c , N ICB /N = 0.34. 

may also be related to the recent specific heat results on 
Ba(Fe 1 _ :c Co :c ) 2 As2^ 

There are two possible contributions to the field depen- 
dence of Tj -1 arising from quasiparticlc excitations in the 
mixed state: localized quasiparticlcs inside vortex cores, 
and delocalized quasiparticlcs originating from nodes in 



the SC gap. Since the applied fields are much smaller 
than H C 2 ~ 52 T^ localized quasiparticles inside vortex 
cores are unlikely to contribute to the relaxation at these 
fields (H/H C 2 < 0.22). In fact, the thermal conductiv- 
ity measurements also show a sizable field dependence 
of kq(H)/T at low fields, attaining nearly 70% of the 
normal state value even at 0.2 H C 2, indicative of the pre- 
dominance of delocalized quasiparticles^ These results 
suggest that the field dependence of T^ 1 is attributable 
to the contribution from delocalized quasiparticles origi- 
nating from nodes in the gap, compatible with the exis- 
tence of a residual DOS in the SC state. 

The existence of a residual DOS in the SC state 
of BaFe2(Aso.67Po.33)2 differs strikingly from other 
iron arsenides, including 1111—^ and 122 superconduc- 
tors f ^ 12 ' 28 instead resembling the cuprates and heavy 
fermion superconductors, which have nodes in their SC 
gap. As shown in Fig. 5, 57 Fe T^ 1 in (Bao.6Ko.4)Fe2As2 
near optimal doping exhibits a T 5 -like dependence well 
below T c r^ The fairly large exponent of ~ 5 in T^ 1 im- 
plies a fully gapped state in (Bao.6Ko.4)Fe2As2, and was 
analyzed in terms of s±-wave superconductivity with two 
gapsii 2 - In the s± models with unitary scatterings ; 29 ' 30 in- 
gap states at zero energy can be induced without severely 
reducing T c and can in principle account for the T^ 1 oc T 
behavior. However, the penetration depth AA(T) oc T at 
low temperatures observed in BaFe2(Aso.67Po.33)2 is in- 
compatible with these models ; 31 ' 32 but is consistent with 
line nodes in the gap^ In a symmetry-imposed nodal 
superconductor such as described by a line-node model 
A = A sin(2(/)), T^ 1 oc T 3 in the clean limit, but im- 
purity scattering induced near gap nodes gives rise to an 
enhanced T^ 1 that is proportional to temperature at low 
temperatures M- P dopants in BaFc 2 (Aso.67Po.33)2 can 
provide a source for nonmagnetic scatterers, consistent 
with the existence of Fermi-level DOS. In fact, the 
data of BaFe2(Aso.67Po.33)2 are well fit by the line-node 
model with 2A /fc B T c = 6 and N Tea /N = 0.34, where 
N les /No is the ratio of a residual DOS in the SC state to 
that in the normal state (see Fig. 5). Therefore, NMR, 
penetration-depth, and thermal conductivity results on 
BaFe2(Aso.67Po.33)2 are consistent with the existence of 
line nodes in the SC gap. 

In the case of nodes imposed by symmetry, a residual 
DOS at zero energy is a natural consequence of the impu- 
rity scattering that broadens the nodes. However, this is 
not clear for the case of "accidental" nodes not imposed 
by symmetry, as in the extended s-wave case because im- 
purity scattering will lift the nodes and produce a more 
isotropic gap^ Our experimental findings impose signif- 
icant constraints on candidate theoretical descriptions of 
iron pnictide superconductivity. 

The strikingly different behavior of Tj -1 in the SC state 
between BaFe 2 (As .67Po.33)2 and (Ba .6Ko.4)Fe 2 As2 is 
surprising because they have comparable T c values and 
similar Fermi surfaces. This implies that high-T c super- 
conductivities with and without nodes are nearly degen- 
erate and that subtle differences can significantly mod- 
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ify the nature of superconductivity ! 18 ' 20 ' 3 — Changes in 
the pnictogen height may lead to nodal superconduc- 
tivity, but its T c should be much lower than that of 
fully-gapped s± states^ However, this is not the case 
in BaFe2(Aso.67Po.33)2- One notable difference inferred 
from band calculation s) 22 ' 25 is that BaFe2(Aso.67Po.33)2 
may be more three-dimensional than (Bao.6Ko.4)Fe2As2, 
suggesting that a three-dimensional model is required. 
Reconciling their comparable T c values and fundamen- 
tally different SC gap functions needs further theoretical 
and experimental works. To clarify this problem, identi- 
fying the nodes' location on Fermi surfaces will be crucial. 

In conclusion, our NMR results show the existence of a 
residual DOS at zero energy in BaFe2(Aso.67Po.33)2- This 



is consistent with the presence of nodes in the SC gap in- 
ferred from a field dependence of T-j -1 , and from thermal 
conductivity and penetration depth results. Our results 
suggest that strikingly different SC gaps appear in iron 
pnictides despite comparable high- T c values, analogous 
phase diagrams, and similar Fermi surfaces. 
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